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Abstract 

The Bd — > ^Ks process is without tree level decay amplitude . Then the SUSY 
contributions may be significant . We calculate the neutralino and gluino penguin 
contributions to the CP violating phase in SU{5) grand unified theory with right- 
handed neutrino SU{5)rn ■ We find that the SUSY contributions change the phase in 

1 Introduction 

In the minimal Standard Model(SM) lepton flavor is conserved . Then Lepton Flavor Viola- 
tion (LFV) and flavor changing neutral current (FCNC) is forbidden at tree level . Therefore, 
the flavor mixing is controlled by the Cabibbo Kobayashi Maskawa (CKM) matrix element 
. So the phase in the CKM matrix is the only source of the CP violation in the Standard 
Model (neglecting the strong CP problem) . 

However, in supersymmetric (SUSY) models there is a new source of flavor mixing in 
the mass matrices of SUSY partners for leptons and quarks . In the context of minimal 
Supergravity (mSUGRA), where SUSY breaking parameters are assumed to be flavor blind 
at the Planck scale, the flavor mixing is controlled by the CKM matrix elements . As a 
result, the maximal deviation from the SM in the CP violating parameters is of the order 
of 10%, while the SUSY contribution to b — > 57 process can be much more important . 
In the GUT scenario, there are additional contributions to FCNC/LFV process from GUT 
interaction^] . As for LFV processe the ^ — > 67 branching ratio is closed to the current 
experimental bound, especially for 50(10) model . Experimental evidences of neutrino 
oscillation indicate the existence of small neutrino mass and large flavor mixing in the lep- 
ton sector . A natural explanation for small neutrino masses is the seesaw mechanismjl] 
. In this mechanism, heavy right-handed neutrinos are introduced . These neutrinos have 
Majorana mass term and new Yukawa interactions . In SUSY models with a neutrino mass 
generation mechanism of seesaw type, Yukawa coupling constants among the Higgs doublet, 
lepton doublets and right-handed neutrinos could induce large flavor mixing effects in the 
slepton sector . In the Supersymmetric Grand Unifled theory (GUT) with right-handed neu- 
trino, above the GUT scale, the right-handed down type squarks couple to the right-handed 
neutrinos . Then, due to renormalization group effect, flavor violation in the slepton sector 
may be transferred to the right-handed down type squark mass matrix^], which affects the 
CP violation in the B decay . 
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The Bd — > process is without tree level decay amplitude . Then the SUSY contri- 

butions may be significant . The SUSY contributions, to the CP violation in this process, 
are gluino, neutralino and chargino . In the Bd — > ^Kg process, gluino and neutralino 
interact with down type squark-quark (fig 10,11 ), whereas the charginos interact with up 
type squarks . So in our model, SU{5)iiN, only gluino and neutralino contributions are 
concerned by fiavor violation in the slepton sector . 

In this context we calculate the neutralino and gluino penguin contributions to the CP 
violation in the process Bd — > ^Ks . Then, the SUSY contributions can be large and the 
phase in the decay amplitude can be as large as 0(0. 1) which is much larger than the SM 
. In our calculation we apply the multimass insertion method ^ . We take into account 
the mixing of the squark eigenstates as well as the mixing in the neutralino sector and we 
consider the case of large mixing between second and third generation, suggested by the 
atmospheric neutrino problem . The calculation of the neutralino contribution is new 
. Previous gluino calculation can be found in the reference equation (A. 16), where mass 
insertion method has been usecfl . In this calculation to obtain the results presented in 
figure (2) by the author, we must take the derivative of the equation (A. 16) of the reference 
P] . Consequently, we have also calculated the gluino contribution, with multimass insertion 
method ( where the left and right squark masses are taken into account) . 

This paper is organized as follows . In section 2 the SU{5) grand unified theory with right- 
handed neutrino SU{5)jiN is introduced . In section 3, we give the explicit expressions of the 
amplitudes associated to the neutralino and gluino contributions to the decay Bd — > ^Kg 
{b — > s$) . Finally, in section 4 we find the numerical results on the SUSY contributions 
to the decay phase for different values on SUSY parameters . In section 5 we conclude . In 
annex, the analytic expressions for the Feynman integrals, which appear in the evaluation 
of the amplitudes, are given . 



2 SU{5) grand unified theory with right-handed 
neutrinos SU{5)rn 

We consider the supersymmetric SU (5)grand unified theory (GUT) with right-handed neu- 
trinos, as an extension of the Minimal Supersymmetric Standard Model with right-handed 
neutrinos (MSSMRN)|S|, 0, • Above, the GUT scale the superpotential is : 

+ TV,; [Yn^j H + ^N, [MN]^J Nj 

where (the quark doublet), $i (the charged lepton doublet) and Ni (the right-handed 
neutrino) are 10, 5 and singlet matter fields of SU{5) in i-th generation, while H and H 
are Higgs fields in 5* and 5 representations, such that : H — {H2,Hc) and H — {Hi, He) . 
Finally, Hi , H2 are the Higgs multiplets of the MSSM and He , He are two different colored 
Higgs multiplets . Here Mjv is the Majorana mass matrix for right-handed neutrinos . For 
simplicity, we chose degenerate Majorana mass matrix for the right-handed neutrinos : 

[M„],j- = Af^„ (5y (2) 

The Yukawa matrix Yu is a complex symmetric matrix, while Yd and Yn are complex . It 
is convenient to work in the basis where Yd is diagonal . Then we have : 

Yu = V^QqYuVq ,Yd = Yd andYN = YnVlOl, 

Yu , Yd andY]^ are real diagonal matrices: 

Yu = diag{yu, Vc, Vt) , Yd = diag{yd, Vs, Vb) , Yn = diag{y^i,y^2, Vvz)- 
^In this method an average squark mass is taken into account only. 



The Q's are diagonal phase matrices : 

Qq = dm5r[exp(i(/)^'^^),exp(i4'^^),exp(i(/)^'^^)] 

such as: 

Moreover, Vq and Vl are unitary mixing matrices parametrized by three mixing angles and 
on CP violating phase . 

At GUT scale (Mqut « 1.9 x lO^^Gey), the GUT gauge symmetry is spontaneously 
broken into the SM one : SU{5)gut ^ SU{3)c x SU{2)l x U{1)y ■ So, at low energy, 

contains (5*(3, 2,^), U (3, 1,— |) and E — (1,1,1), quark and lepton superfields . $ 

includes D (3, 1, i) and = (1,2,— i) quark and lepton superfields . iV is a singlet of 
SU{3)c X SU{2)l X U{1)y ■ The representations for SU(3), SU(2) groups and U{1)y charge 
are indicated in parenthesis . Then the GUT multiplets in terms of MSSM field are: 

^ {Q, V^QqU, OlE}, , ^ (D, e[i},. (3) 
The superpotentialQ becomes: 

W = Q^[V^Yu],jUjH2 + Q^[YDhDJHl + E^IYe^L^Hi 

+N,[YNVLhL,H2 + \M,^N,N, ' ' 

In the SU{5) limit, the naive GUT relation is predicted Ye = Yd at the GUT scale . 
Although this relation gives a reasonable agreement for mi, and rrir , it is well know that the 
mass ratio of down-type quarks and charged leptons in the first and the second generation 
can not be explained in this wajH . Then, in Q, the unitary matrix Vq becomes the CKM 
matrix :Vq ~ Vckm ■ The left-handed neutrino mass matrix induced by seesaw mechanism 
is : 

with = 2(< /ii >2 + < /i2 >2) « (246Gew)2 and tan(/3) = . Then Vl plays the role 

of the neutrino mixing matrix and the neutrino Yukawa coupling is : 

= X f{VL,m,,taniP)). (6) 

For a large M^^^ the equation lO becomes non-perturbative below the Planck scale . In terms 
of MSSM fields, the third term in the superpotentiald can be written as : 

W = NdYNVLQLhD.H, + = ^ y,. [14],^ exp{i^f^)N,D,H, + (7) 

The phases t/)*^'^^ and 0'-'^^ remain in the colored vertices . The effect of the colored Higgs is 
negligible at tree level, but affects the soft SUSY breaking mass parameters through the RG 
effect and generates off-diagonal elements in the mass matrix of down-type squarks . In this 
model we assume the universality of the scalar mass (too) at the reduced Planck scale M, ~ 
2.4 X lO^^Gew . Thus, we adopt the model so-called mSUGRA . Due to the neutrino Yukawa 
matrix, non-vanishing off diagonal elements are generated . The off-diagonal elements in 
the mass matrix of right-handed down-type squarks (D) are approximately : 



[m|.J., ^ - exp[-^(0f ) - cj^f)] vl, [Vnk^[VL]k, x (Sto^ + al) \og{—^) (8) 



^One possibility to remedy this defect is to introduce higher dimensional operators suppressed by the 
Planck scale. 



where ag is the universal A parameter and R stands for right component . For the left- 
handed down- type squarks, the off-diagonal elements in the mass matrix are generated by 
the top Yukawa coupHng, such as: 



[m'~ ^ ~^ynVSKM]u[VcKMh X (3mg + al) x (31og(— ^) + log(-r^)) (9) 

where L stands for left component . These off-diagonal elements mass matrix [m~ ]i3 and 

[m~ ]i3 are coefficients of AB ^ operators, they change the standard model predictions 

to the mixing and the decay of the B-mesons . 

In this context we calculate the neutrahno and gluino penguin contributions to the CP 
violation in the process : Bd — > ^K^- 



3 Neutralino and gluino penguin contributions to 

The operators contributing to Bd — > [H] have a structure like (s6)(ss), then they are 
induced at the one-loop level in the standard model . Therefore the SUSY contributions 
may be significant . The SUSY contribution to the decay amplitude for Bd — > ^Kg comes 
from the penguin and box diagrams . In our analysis, we only take into account the chromo- 
dipole contribution because, as is mentioned in reference|2|, box and color-charge form factor 
contributions have the same size but the opposite sigr|f| . Consequently, the AB — 1 effective 
Lagrangian considered takes the following form: 

Leff = mtCS^^'T^,r[rn'']PLb'G;,, + iL^R) + h.c (10) 

where nib is the bottom quark mass,T^ is the SU{3)c generator, C^^ is the gluon field and 
a and b the color indices. 

We have calculated the neutralino and gluino contributions to the coefficient C^*^ . 
We apply the multimass insertion method .The calculation of the neutralino contribution is 
new . Previous gluino calculation can be found in the reference jH] equation(A.16), where 
mass insertion method has been used . The results are very large and to obtain the results 
presented in figure(2) by the author, we must take the derivative of the equation(A.16) 
of the reference P] . Consequently, we have also calculated the gluino contribution, with 
multimass insertion method, to the chromo-dipole operator . 



3.1 Neutralino contribution 

We assume neutrino and quark mass hierarchy and y^^ ^ Vi^i ■ From the atmospheric 

neutrino observation we take : 

-for the neutrino tau mass :m^^ ~ 0.07ev 

-and large mixing angle :Vl32 = 

Thus, from the equations ijSl and lO we have : 



-iexp(-*(4^' -4^)))y2jF/]33[l^L]32 X i3ml + al)logi^^) (11) 

and 



[m| ]32 ^ -^y^VSKMUVcKMhs X (Sm^ + ag) x (3 log(-^) -t- log(|^)). (12) 
''^ 8tt^ Mgut Mweak 



^In the neutralino case, the wino and higgsino couplings have the same sign as gluino ones . Then 
neutralino and gluino contributions have the same behaviour . For chargino it is diflferent . 



These two terms contribute to the decay amphtude: 



MifJ^^o ^K'\Leff\B, > . (13) 

The right-right non-diagonal termfentlTll depends on the neutrino phase and is dominant 
compare to the left-left one (eqEl because |[VcA:M]ts| ^ |[Vl]32| • So, we only consider 

( SU SY) 

the right-right contribution to the Mg^_^^-^a decay amplitude . Feynman rules used for the 
calculation are summarized in figures(10)-(ll) . In annex, we have the Feynman integrals . 
In the following we note:(m^)32 = [to— ]32 .The three kind of diagrams contributing to the 

decay are illustrated in figure 1 and figure 2 . 

From the diagrams illustrated in figurel we obtain : 
+for the diagrams 1(a) and 1(b) 



rp V2g2g3 tan(6^)x (mi,+ms)Mi /i (m )32 

^'^a--b — 3xl6Tr2[m2, -mi ][Mf-/i2] 

bR sR 



^ J2n=0 



FCi{X~ )+FCi{X~ ) FNi{X~ )-FNi{X~ ) 



(14) 



2 J "V for n = 2 

with M„ = < 9 " . and X~ — =2- and X~ — , where Mi , a are the bino 

" ' for n ~ 1 M„ A' M„ ' ' ^ 



and higgsino masses, mi and mA stand for the diagonal components of the bottom and 

bR sr 

strange squark mass matrix (m^ is the strange quark mass) . The coupling expressions are : 
52 = 4"''ri^^9'^) ' 93 ~ ^^4T^' 'I outgoing gluon momentum and p the bottom 

quark one . 

+on the diagrams l(c)-l(d), there is a flip of chirality on the bottom leg then. 



Ic- 



gjga TOfc (to )32 

^ hR SR^ 



TOsTOf, 



2M2cos2(/3) 9 



tan 



n=0 J''^ n 



(15) 

The functions FC\{X) , FC2{X) , and FNi{X) are given in annex. 
The contributions, illustrated in flgurel(e-f-g-h), are proportional to TOf, due to the left- 
right diagonal insertion : 



-aoTOoTOb — TOh/itan(/3). 



(16) 



Then from these diagrams, we get: 
+for the diagrams 1(e) and 1(f) 



rp _ 9293 "»I.flb ("^ )32 



bR bL 



Ml tan 



9 



2AfJ cos^(/3) 



y^^ _(-ir X J: 



2xMiXFCi{X~ ) 2x/j^xP'Ci(X~ ) 

=3 ^ -h =2 - 3 X FCsiX- ) - 3 X FCsiX- ) 



(17) 



where 



„ f mi for n — 

mi- f or n — 0,1 , ■j'''^ 

?7ii forn = 2,3 | ^« 

f>i ?Ti4 for n — 3 



(18) 



+for the diagrams 1(g) and 1(h) 



^3 (-1)" „ 1 



2(^^+2^Afi) 



(19) 



2(M^+2pM: 



with: 



'■n2 



mi /or n = 

r7i4 /orn = 0. 1 , 

nii- for n = 2, 3 



for n — ?> 

bR ■' 



(20) 



The functions FC3{X),FC4{X),FC5{X) and FC6(X) are written in annex . 
Finally, the contributions drawn on figure 2(a-b) and 2(c-d) are identical with the last 
ones (fig.l:e-h), but they are proportional to due to the left-right diagonal insertion : 



7^2 



'■LRs 



— —aomoTris — m^/i tan(/3). 



(21) 



Therefore, the expressions T2a~b and T^c-d are obtained by interchanging : ^ , hs. ^ 
Si? , 'sr ^ sl and mi, ^ iris mj^jn, ^ tt^^rs in the equations IjlTL 1181 1191 12011 . 

3.2 Gluino contribution. 

Due to gluino-gluon interaction and the multimass insertion method, the eight diagrams 
(four in the mass insertion method case ) contributing to the chromo-dipole are illustrated 
in figure 3 . Due to the fiip of chirality on the strange leg, the amplitudes associated to the 
diagrams(3c) and (3g) are proportional to rus ■ Then, these two contributions are negligible 
comparatively to the others . 

From the diagrams illustrated in figureS we obtain : 



[3a 



ri=0 *J n 



2N, X 167r2[mi -mi 



Mi 



with Xg = , M~ is the gluino mass, Nc — 3 and 

rri~ for n ~ 
Ml = { rf^^forn = 1,2 
for n — "i 

bL 



(22) 



(23) 



and 



_ 93M~{m )32'ni.flfc 3 (-1)" 

^ 3fc - 2AfcXl67r2pSi ~n& ] ^ ^n=0 mi -nit 



M„ 



2 Mi 



M, 



bR bL 



■'-FC^{Xg)~iFC3{Xg) 



(24) 



^'1=0 —rr, 



3d 2Afe X 167r2 [rfCd, -mi. ] 

sR bR 

2 X Mi. 



Ml 



^FC,{Xg,)-3FCsiXg,) 



(25) 



with 



nii for n = Q 



'n3 



for n — OA 

TnA for n = 2, 3 
bR ■' ' 



si?. 



M 



n3 



From the others diagrams we deduce : 



T, 



3e 



2 X 167r2[mi -mA 1 

^ bR. sR' 



E 



nrd- for n — 1,2 
TT^E for n = 3 



(-1)" 



M„ 
3 



3/ 2xl67r2[mi -mi 1 Z^ri=0 m^,,-mi 



^Gi(X,2) + 3FC3(X32) 



'3/1 



7^ 



(-1)" 



m„ , — mi 



=^Gi(Xg4) + 3G2(Xg4) 



with 



TTir^ for n = 0, 1 
mA for n = 2, 3 



tttE for n — Q 
_J>R 

m~ for n — 1,2 

si 

mi for n = 3 
si? •' 



(26) 



(27) 



(28) 



(29) 



(30) 



Due to the multimass insertion method, all the amplitudes depend on the squark left 
and right masses (contrary to the mass insertion method case where an average squark mass 
contribute ) . Now, we can build the SUSY contribution to the decay amplitude : M^^^^j^g- 



4 Numerical results on SUSY contribution to the decay 
phase 

In the decay amplitude, we take into account the chromo-dipole term only . Then, with the 
factorization approximation pi] -^2], the decay amnlitude lfT^jl takes the following form : 

Mb,-,4>ko = 2 X (Pse^)m2/^Ff ^1(1 - ^) [gsKoMCg'' + {L ^ R)] (31) 



where < ^(P^, e0)|sVslO mj/^e'' 



and < K^{PK)\T-i''h''\Bd{PB) >= Ff^iPg + PrY + F^^{Pb - Pk) 



(32) 



The Kdm coefficient is induced from the hadronization of the chromo-dipole moment 
From the quark model and the heavy quark effective theory, we have : 



KnM ~ ^ 



^ + 0{ml/ml) , 
where q is the momentum transfer in the gluon line such that : 



9 



2 1. 

1 = o("l 



1 2 



(33) 



(34) 



and Kdm ~ l-'i for as{mz) — 0.118 , aem(™z) = 1/128 and — 5Gev . In the decay 
process Bd — > the CP violation is determined by : 



/total 



imix I _ 



M. 



{SUSY) 



So, the SUSY contribution to the decay phase is expressed by : 



(35) 



AfSUSY 



Bs^<t>KO 



(36) 



where is the maximal possible correction to the decay phase such that: 



4^) 



t>2 — arg 



Mi" 



Bd^<j>KO 



TT 

2' 



(37) 



The standard model contribution ^21 has the following form : 



MSM ^ _^VtbV,lm,lUe^ x 2P^F^^ 
C+^g.^l^(m, + 4("l + ^)-™,(l-2 ^ 



X 



4m 



,jnb\ 



(38) 



with C = C3 + C4 + C5 + ]^ (C3 + C4 + Ce) - i Cy + Cg + Clo + ^{cs + C9 + Clo) 



1 



and the 



are the Wilson coefHcients . The equation (l36)l is very sensitive to the standard model contri- 
bution . So we choose to take into account the next leading order QCD corrections, with the 
renormalization scale :/z = mi, = 5Gev,aemi'mz) = 1/128, as{mz) = O.llSandmt = ll&Gev 
. So, the Wilson coefficients have the following values : 



C3 = 0.0174 
C7 = 1.398 X IO- 



CS 



C4 -0.0373 
= 3.919 X 10-4 



C5 



0.0104 
-0.0103 



cio 



C6 = -0.0459 
= 1.987 X 10-3, 



and the DPO coefficient at the two loop level is cn — —0.299. 
Concerning the SUSY contributions, we have considered: 
+the neutralino contribution 



Bi 



= 2{PBS^)mlUF 



BK 



x^(l-4)53i^. 



DM 



with Tx = Tia-b + Tic-d + Tie-f + Tig-h 



T2a-b + T2c-d and 



(39) 



M2i 

Ba^4,K'> 



MSM 

Bs^<l>K« 



-the gluino contribution 

^„ = 2{PBe^)mlUFf''^ 



Bd 



with T~ = Taa 



+ T3b 



L3d 



3e 



T'- 



3/ 



T3/1 and 



7V2 



)g3KDM X X Tp 



B,->4,K0 



(40) 



(41) 



(42) 



For each case and different Kdm values, we have calculated the decay phase . From the 
atmospheric result, we take the tau neutrino mass : m^^ — 0.07ey and Vl32 = ""Tf ' '^^^ 
right-handed neutrino mass is constrained by the neutrino Yukawa coupling eq® and the 
Br {fj, 67) experimental bound[2] ■ So we choose Mi,^ — 2 x 10^'^GeT^ . In our numerical 
analysis we have used the suspect^! program by which is possible to obtain the sparticle 
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Too 


200 


300 


400 


500 


600 


700 


800 


1000 


1200 


1500 


Ml 


81.55 


81.79 


82 


82.27 


82.56 


82.81 


83 


83.48 


83.9 


84.3 




309.1 


327.7 


352.5 


381.5 


413.7 


448.2 


483.5 


564.5 


640.9 


761.1 


TO— 

SR 

TO— 

ST, 


466.7 
484.6 


515.5 
531.7 


577.6 
592 


648.8 
661.7 


727.5 
738.9 


810.5 
820.8 


897.2 
906.5 


1078 
885.6 


1264 
1271 


1550 
1556 


m~ 

bR 

mr 

br. 


466.4 
433.6 


515.1 
466.6 


577.2 
510 


648.3 
561.1 


726.9 
619.1 


809.9 
681.3 


896.5 
747.4 


1077 
885.6 


1263 
1032 


1549 
1259 


M~ 

Li 


511.7 


518.2 


525.1 


531.2 


538.4 


544.6 


550.6 


561 


570.6 


582.6 



Table 1: Sparticle mass spectrum (en GeV) for /i > 0,toi/2 ~ 2QQGeV and tan/3 = 3 



Too 


200 


300 


400 


500 


600 


700 


800 


1000 


1200 


1500 


Ml 


81.48 


81.75 


82.06 


82.32 


82.48 


82.74 


82.96 


83.4 


83.75 


84.31 




263.7 


265.6 


267.7 


270.3 


272 


274.4 


277.3 


284.9 


278.2 


265.5 


m— 

sr 

TO— 

sr, 


465.8 
483.7 


514.9 
531.1 


577.1 
591.6 


648.7 
661.6 


726.7 
738.1 


809.9 
820.2 


896.5 
905.8 


1077 
1085 


1264 
1270 


1550 
1555 


TOr 

bR 
br, 


462.9 

435 


511.5 
469.5 


573 
514.4 


643.8 
567.2 


721 
626.1 


803.4 
690 


889.2 
757.3 


1068 
899.1 


1253 
1049 


1537 
1280 


M~ 

Li 


511.3 


518.3 


525.3 


532.2 


538.1 


544.4 


550 


560.8 


570.2 


583.3 



Table 2: Sparticle mass spectrum (en GeV) for ^ > 0,TOi/2 = 200GeT^ and tan (3 — 10 



mass spectrum for different values on SUSY parameters that is :too , TOi/2 , tan/3, sign{fi) . 
In this program RG evolution of parameters is taken into account . To fix the gluino mass 
we take toi/2 — 200 Gev and sign{fi) positive . These values are given in tablesEl[2l|21and 
are consistent with Br (/i ej) experimental bound . 



We plot in figures 4-6, the neutrahno contribution to tan(A(/)^^^^^^);^ = 

and gluino ones ea lj42|l on figures 7-9 . 

From the results obtained, we remark the following . 



M2i 



(1) Because of our use of multimass insertion method, the results depend on the gluino, 
wino, higgsino and squark left and right masses (tables [2 El ) • 

(2) The gluino and neutralino contributions have the same behaviour and the upper values 
are obtained for too = 500 Gev (figures 4-9) . 

(3) As is shown on fig. 4 and fig. 7, these two contributions are very sensitive to tan/3 and 
Kdm values . They increase with the enhancement of these two parameters . 

(4) The dependency on of the decay phase is presented on fig. 6 and fig. 9 . The phase 
becomes more important with increasing ag parameter . 



Too 


200 


300 


400 


500 


600 


700 


800 


1000 


1200 


1500 


Ml 


81.53 


81.8 


82.02 


82.27 


82.55 


82.81 


82.99 


83.5 


83.81 


84.24 




258.8 


258.8 


258.7 


258.7 


257.8 


256.2 


255.1 


253.4 


237.7 


190.8 


m~ 

SR 

TO— 


466 


515.2 


576.9 


648.5 


727 


810.1 


896.7 


1077 


1264 


1550 


483.9 


531.4 


591.4 


661.3 


738.5 


820.4 


906 


1085 


1270 


1555 


mr 

bR 

mr 

br. 


444.1 


488.4 


544.3 


609.2 


680.7 


756.3 


835.5 


1001 


1172 


1436 


425.8 


457.8 


499.2 


548.6 


604.3 


664 


727.3 


861 


1002 


1222 


M~ 

Li 


511 


517.9 


524.2 


531 


537.9 


544 


549.6 


561 


569.8 


582.1 



Table 3: Sparticle mass spectrum (en GeV) for n > 0,toi/2 = 200GeT^ and tan/3 = 30 



(5) In the gluino case, thanks to the gluino mass (500 GeV), the decay phase is significant 
for large tan/? and moderate K^m values . Indeed, from fig. 8, for oq = 0, tan/? = 30 and 
Kdm — 1-5 the upper value of the decay phase is 4.10"^ . This value is greater than the 
uncertainties in the SM calculation O(10~^) |S] ■ Moreover, as is shown on fig. 9, for oo = 10 
the gluino contribution is more competitive with the SM one because the upper value is 0.1 . 
Our results are consistent with those of the reference[9J (figure 2) . However, our results are 
much smaller due to multimass insertion method and the value of the Afj,^ = 2 x \Q^'^GeV 
(instead of M^^ = 5 x lO^^GeV in the reference[ni) . 

(6) For the neutralino case, due to the small values of n and Mi (tables 12 El ) , the decay 
phase (fig. 4, 5) cannot be competitive with the gluino one for ao less than 250 (even for 
large tan/? and Kom values) . The highest value (10"'^) is very small in comparison with 
the uncertainties in the SM model calculation . Nevertheless, for large ao (ao > 250) and 
tan/? values, we obtain results (fig. 6) greater than 10^^ . Then, the neutralino contribution 
changes the phase in the process Bd — > . 



5 Conclusion 

In the context of SU{b) grand unified theory with right-handed neutrino, we have calculated 
the neutrahno and gluino penguin contributions to the phase in Bd — > ^Ks ■ In our 
calculation, we apply the multimass insertion method . For convenient choices of SUSY 
parameters, the neutralino contribution to the decay phase is greater than O(10~^) . For 
the gluino contribution the decay phase is the order of 0.1 . Then, the SUSY contributions 
change the decay phase in the process Bd — > ^Kg ■ 



6 Annex 

Now, we give analytic expressions of the Feynman integrals which appear in the evaluation 
of the amplitudes : 



Q^(^X) — ~4-llX-14X log(X) + 16X^-4X^ \og{X)-X-' 



^2 ^ ^2 1 

2X(X~l)i- 



G2{X) 



-38-24 log(X)+81X+18X log(X)-54X'' + llJi:^ 
36(X-l)'i 



i*05^Aj - 3(X-1)B 

T?n (V\ — -1+6X-3X^+6X^ \og(X)-X'^ , 

TPn (Y\ 5+2 log(X)-4X+4X \o^(X)-X^ 
tCi{X) - 2(X-l)i 

FC2{X) = ~l+~-^-3X^+6J^^ log(JC)-2X^ 



6(Jf-l)* 



_ -2-3X-6Xlog(X)+6X^-X ^ 
~ 6(X-1)* 



_ l-8X-12jy^ log(Ji:)+8X^-X-* 
— 12(X-1) 



_ l+2Xlog(X)-X" 
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Figurel:Neutralino contributions with gluon (g) and squarks (s^, &) coupling to right-right 
non-diagonal mass insertion (m^)32, with left-right diagonal mass insertion m^jf^jj, Mcos = 
\pim^ cos(/3) and Mg-m = V^fn^, sin(/3) . 5° and are bino. higgsino components of the 
neutralino (L stands for left and R for right) . Mi and // are the bino and higgsino masses 
respectively. 




Figure2:Neutralino contributions with giuon (g) and squarks {s,b) coupling to right-right non- 
diagonal mass insertion (m^)32, with left-right diagonal mass insertion Mcos = V^'m^ cos(/3) 
and Msin = V^rriy, sin(/3) . S° and are bino, higgsino components of the neutralino (L stands 
for left and R for right) . Mi and fi are the bino and higgsino masses respectively. 
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Figures :Gluino contributions with gluon [g] and squarks coupling to right-right non-diagonal 



mass insertion . 
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figure 4: 10-^tan(A$^"^^_^^Jx(^o) forKDM=l, M^^ = 2W^GeV, mi/2 = 
200GeV and qq = 
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figure 9:tan(A$^"^^_^^Jg(mo) for tan(^) = 30, M^^ = 210^^GeV, mi/2 = 
200GeV and KDM''=1.5 
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FigurelO:Gluonic Feynman rules, where T^^ is the SU{3)c generator, g the squark d,^ the gluino,g the 
gluon and Pl and Pr are the chirahty projector. 
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Figurell:Neutralino Feynman rules. is the squark down.i7° and stand for higgsino and bino 
component of the neutrahno-Pj, and Pr are the chirality projectors. 



